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Abstract 

We investigate electron charge decoherence in a laterally-coupled single-electron semiconduc- 
tor double quantum dot through electron-phonon interaction. We analytically and numerically 
evaluate the relaxation and dephasing rates due to electron coupling to both acoustic and optical 
phonons, and explore the system parameter space in terms of interdot distance, strength of single- 
dot confinement, and inter-dot coupling strength. Our numerical results show that the electron 
scattering rates are strongly dependent on the strength of the electron confinement and the size 
of the system. In addition, although the most dominant factor that determines the charge deco- 
herence rate is the energy splitting between the charge qubit states, the details of the double dot 
configuration is also very important. 

PACS numbers: 03.67.Lx, 73. 21. La, 85.35.Be, 63.20.Kr 



I. INTRODUCTION 



Ever since Peter Shor showed that a purely quantum mechanical computer can be used 
to achieve exponential speedup (compared to classical computers) in solving the prime fac- 
toring problem, 1 there has been widespread interests in the study of quantum information 
science in general, and in building a practical quantum computer in particular.^* 3 - Some of 
the most prominent proposals are based on solid state structures, whether superconducting 
nanocircuit a 4 i 5 i 6 i 7 or semiconductor nanostructures such as quantum dots and regular donor 
arrays . 8 ! 9 ! 10 ! 11 ! 12 ! 13 A major presumed advantage of these systems, especially the semicon- 
ductor artificial structures, are their potential scalability, supported by the powerful and 
advanced semiconductor industry. 

The great interest in semiconductor-based quantum computer architectures has prompted 
extensive studies of a variety of physical properties of nanostructures that are relevant to 
single electron spin or charge degree of freedom.— In the present paper we focus on the 
decoherence properties of a charge qubit in semiconductors . 10 1 12 > 13 For a charge qubit based 
on a single electron in a semiconductor double quantum dot (QD), charge decoherence 
has two important channels: Coulomb interaction to the background charge fluctuation 
and electron-phonon interaction . 15 i 16 i 17 i 18 i 19 i 20 i 21 The former decoherence channel is widely 
present in all nanostructures. Presumably it is due to coupling with defects in the system 
and is thus extrinsic. The latter, however, is intrinsic to any solid state host material for 
the charge qubit. 

In the following we present our study of charge qubit decoherence caused by electron- 
phonon interaction in a horizontally coupled two-dimensional GaAs double quantum dot. 
In Section II, we identify the electronic states we are interested in, and clarify the relevant 
types of phonons and electron-phonon interactions involved. We then derive the various 
relaxation and dephasing rates for double-quantum-dot-trapped single electrons. In Section 
III and IV, we show our results on charge relaxation and dephasing rates for a variety of 
configurations and states, and discuss the physical pictures and implications. Section V 
presents a summary of our results and our conclusions. 
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II. THEORETICAL DESCRIPTION OF ELECTRON-PHONON INTERACTION 
A. A single confined electron in two coupled quantum dots 

In this study we consider gated lateral quantum dots in a AlGaAs/GaAs/AlGaAs quan- 
tum well (QW). The growth direction (z direction, or vertical direction) confinement is 
due to the higher bandgap of the barrier material of AlGaAs. The lateral confinement is 
produced by the electrostatic potential from surface metallic gates. In general, the vertical 
direction confinement length (~ 10 nm) is much smaller than the lateral confinement length 
(50 nm), so that we can safely treat the dynamics along vertical and horizontal directions 
as decoupled. The system Hamiltonian (within the effective mass and envelope function 
approximations) is thus 

n = n ll+ n z ^ (i) 

where the growth direction component Ti, z takes the form 

n z = -h z ^—d z + v e (\z\ - l z ) . (2) 

2 m*(z) 

Here m*(z) is the electron effective mass and Vq the offset between the band edges of the 
GaAs well and the AlGaAs barrier. For simplicity, we take the ^-direction wavefunction as 
the wavefunction of an infinite QW (Vo — > oo). In this work, we do not consider excitations 
along z-direction because of the much higher excitation energy (compared to the lateral 
direction), so that the ^-direction wavefunction is always given by ip z (z) = Acos(irz/2L z ) 
where A is a coefficient to be determined by normalization and 2L Z is the width of the QW. 

The lateral confinement is assumed to be parabolic for a single QD, so that a single 
electron Hamiltonian in the lateral direction is 

where ujq describes the strength of the harmonic confinement in the x — y plane. The total 
electron wavefunction can now be written as a product of 

V>(r) = (r||) V* (*) (4) 

In the case of a single QD, the two-dimensional (2D) one-electron wavefunctions are essen- 
tially 2D harmonic oscillator function o 22 ' 23 and are described in terms of the principal quan- 
tum number n — 0, 1, 2, ... and the angular momentum quantum number m — 0, ±1, ±2, ... 
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as 



7r/ 2 (n + | m |)! 

where (p 2 ) are the Laguerre polynomials, and p = |rn|// is a scaled radius, with I 



Ti/m*ojQ. The corresponding eigenvalues are 

E nm = (2n + \m\ + I) huj . (6) 

For two QDs that are horizontally coupled, we use a simple in-plane confinement of two 
parabolic wells separated by an inter-dot distance 2a: 

V c = ^m*ul min{(a; - a) 2 + y 2 , (x + a) 2 + y 2 } (7) 

The single electron wavefunction for the lateral direction is, in general, given by a superpo- 
sition of the single-dot wavefunctions: 



*ll) = £ff* ^u) + D k Vfi*), (8) 

k 

and the total wavefunction of the system of the coupled QDs is 

tf(r) = *|| (r,|) ^ (z) (9) 

Notice that for charge qubits there is only a single electron in a double dot, in contrast with 
spin qubits, where each quantum dot has an electron and double dot is only for two-qubit 
operationsj&24 In the present study, the wavefunctions for the coupled-QD are calculated 
numerically by direct diagonalization, using reasonable parameters of a GaAs QW. 

B. Electron- Acoustic phonons coupling 

In a polar semiconductor like GaAs, electrons couple to all types of phonons. More 
specifically, in GaAs electrons couple to longitudinal acoustic phonons through a deforma- 
tion potential, to longitudinal and transverse acoustic phonons through piezoelectric interac- 
tion, and to the optical phonons through the polar interaction.— The deformation potential 
electron-phonon interaction is given by 



Hu 1 2 Pm v^j lqb(q)( a «i+ o U). ( 10 ) 
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m/m e 


D (eV) 


c s (m/s) 


p (Kgr/m 3 ) 


ei4 (V/m) 










GaAs 


0.067 


8.6 


3700 


5300 


1.38xl0 9 


12.9 


10.89 


36.25 


33.29 



TABLE I: Material parameters (m e is the electron mass). 



where D is the deformation constant, p m is the mass density of the host material, V is 
the volume of the sample, a q and aL q are phonon annihilation and creation operators, and 
p(q) is the electron density operator. Table I presents the material parameters used in our 
numerical calculation. 

Electrons can interact also with longitudinal and transverse acoustic phonons through 
piezoelectric interaction. This type of interaction is essentially due to the lack of symmetry 
of in the crystal, thus for materials like Si, which has crystal inversion symmetry, the piezo- 
electric interaction is not present. On the other hand, the crystal of GaAs lacks inversion 
symmetry, so piezoelectric interaction is nonvanishing. The electric displacement T> is re- 
lated to the electric field S, strain S, and the permitivity tensor e in a piezoelectric crystal 
by2& 

V i = E e iJ £ j + E e iklSkl (11) 
i k,l 

where the third rank tensor e^i is the piezoelectric constant tensor. As a result, the electron- 
phonon coupling due to the piezoelectric effect is 25 



E 



1/2 



A€(<?)p(q)K 



,t 



(12) 



q ^PmVuJ^J 

where A denotes polarization of the acoustic phonons. In the case of zincblende crystals 
there is only one independent and nonvanishing piezoelectric constant: = e^^ = e^, so 
that the matrix element is given b y 25l27l2S 



M p x z (q) = 2e eu {q x q y (z + QyQzix + q x qzi y 



(13) 



where £ denotes the unit polarization vector and e is the electron charge. 

Notice that the Hamiltonian for deformation potential and piezoelectric interaction are 
real and imaginary, respectively, which allow us to investigate separately these interactions 
and calculate the total contribution by simply adding up the two rates.— 



5 



C. Electron-Optical phonons coupling 

Even though the electronic energy involved in the study of charge qubits is generally 
quite small (a few meV) compared to the optical phonon energy (~ 36 meV in GaAs), 
we will demonstrate later that electron-optical-phonon interaction does play a role in the 
decoherence of electron orbital states. For the purpose of this calculation, we will use 
the simple polar interaction in the bulk, neglecting the more intricate details involving 
heterostructures . 29 i 30 i 31 The electron-phonon interaction due to LO phonons is thus given 



in which ujlo is the longitudinal optical frequencies, e s and are the static and high 
frequency dielectric constant. 

Now that we have identified all the electron-phonon interaction involved in our system, 
we are ready to study what concrete forms they take in a double quantum dot system, and 
evaluate the dephasing and relaxation rates for an electron in a double quantum dot due to 
its interaction with the phonons. 

III. RELAXATION AND DEPHASING RATES DUE TO ELECTRON-PHONON 
INTERACTION 

A. Electron-phonon coupling in a double quantum dot 

Before performing numerical evaluations of the electron charge relaxation and dephasing 
rates, we first clarify the physical picture of charge decoherence in a double quantum dot. 
As an example and without loss of generality, let us examine the deformation potential 
electron-phonon interaction given by Eq. in which p(q) is the Fourier transform of the 
electron density operator: 



by25 




(14) 



and 




(15) 




(16) 
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where k and 77 are indices of electronic states, c K and cj, are electronic annihilation and 
creation operators for the /c-state, while <f> are the electron wavef unctions. In the context of 
a coupled double QD, we can choose the double dot eigenstates as the basis for the single 
electron. For the two lowest-energy double-dot states <3>±, which are chosen as the charge 
qubit basis states, k and 77 take the values of + and — (from now on, these two states will 
be equivalently referred to as the ground and first excited states of the double dot, or the 
charge qubit states). The electron-phonon coupling Hamiltonian can then be conveniently 
written in this quasi-two-level basis in terms of the Pauli spin matrices a x and a z (where 
spin up and down states refer to the two electronic eigenstates): 

( ft V /2 

q V^Pm^Wq/ 

A-(q) = He iq - r |+>, 

M<0 = ^((+|e Jq - r |+)-(-|e iq - r |-)) • (17) 

Since the basis for the quasi-two-level system are electron eigenstates for the double dot, 
the term proportional to a x above leads to transition between the two electronic eigenstates 
and causes relaxation. On the other hand, the term proportional to a z does not mix the 
electronic states, so that it only causes fluctuations in the energy splitting between the two 
electronic levels. Therefore it only leads to pure dephasing between the two electronic charge 
states, but not to relaxation. 

For a discussion of the qualitative behavior of the electron-phonon coupling in a double 
dot, we first analyze the simple situation where the two dots are well separated and not 
strongly biased, so that only the two single-dot ground orbital states are involved. The 
relevant single-electron double-dot states are then 

$+ = a0 A (r) + b<f) B (r) ; $_ = b<p A {j) - a0 B (r), (18) 

with 4>A(B)( r ) = f( r — R-a(b)) ^o( r ), where (p(r) is a slowly varying envelope function, and 
the Bloch function at the conduction band minimum (k = at T point) is equal to the 
periodic part Uq(t). Though we have chosen the envelopes if centered at each well to be 
identical, they could as well be different, as is generally the case for quantum dots. For 
small energy splittings between the $± states, the fast oscillatory Bloch function Wo( r ) can 
be integrated separately, so that the matrix element A r can be written as^: 

A r {q) = (ab* - a*6e iq R ) J dr e^Xr)] 2 
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+ (H 2 - M 2 ) / dr e iq >(r)^(r - R) . (19) 

Here the first integral is an on-site contribution modified by the phase difference e iq R be- 
tween the two dots, while the second integral is a two-dot contribution that is generally 
much smaller because of the small interdot overlap. 

The dephasing matrix element A v can be similarly calculated and the result ist^ 

^(q) = z (|&| 2 -|a| 2 )e iq - R / 2 sin^ J dr e^ r [ V (r)] 2 

+ {a*b + ab*) J dr e iq 'V(r)^(r - R) . (20) 

Here the prefactors \b\ 2 — \a\ 2 and a*b + ab* are for intradot and interdot integrals, exactly 
the opposite to those in Eq. (fT9*j) . 

Equations and EO] clearly demonstrates that electron-phonon interaction induced elec- 
tron charge decoherence is dominated by relaxation when \b\ ~ \a\ so that A v is small , 2(|>i21 
and by pure dephasing when \b\ and \a\ are very different (so that, for example, ~ 1 
and \a\ ~ 0). 21 In other words, if the energy levels of the two quantum dots are close to 
resonance, relaxation matrix element is much larger than pure dephasing matrix element; 
while when the two single-dot levels are biased, relaxation is suppressed. These qualitative 
trends will persist even when we consider more realistic electron eigenstates as studied in 
section II. A. The argument here is based on the assumption that the off-site contributions 
are small because of a small wavefunction overlap between the two dots. Obviously, if the 
overlap is larger, the above trend becomes weaker. 

B. Relaxation rates 

The electron relaxation rates associated with phonon emission (or absorption) can be 
evaluated using Fermi's golden rule: 

T = ^^\(¥ F \r)\H int \¥ I \r))\ 2 6{E F -E I ±E (l ) 
n q 

(N B (E (l ,T lat )+±±^ (21) 

Here labels T and 'F' refer to the initial and final electron orbital states respectively, the 
plus (minus) sign denotes emission (absorption) of a phonon, and Nb is the Bose-Einstein 



S 



distribution for phonons with lattice temperature T\ at (our calculations presented in this 
paper are all done at T\ at = 0, when phonon absorption can be neglected.). 

Due to their large energies (~ 36 meV for LO phonons in GaAs), optical phonons do 
not contribute to electron orbital relaxation in a quantum dot except for the highly excited 
states, which are irrelevant for the context of charge-based quantum computing. Therefore, 
the only phonons that contribute to electron relaxation here are acoustic phonons.— 



C. Pure dephasing rates 

Relaxation is not the only way charge qubits can be decohered. If the energy difference 
between the two charge states fluctuates, phase information will get lost and decoherence 
occurs. Such pure dephasing (in the sense that no transition occurs between the two charge 
states) due to a bosonic bath has been calculated before . ?9i?7 Pure dephasing due to electron- 
acoustic-phonon interaction has also been evaluated for spherical quantum dots and donors 
in semiconductors.— The density operator of an electron in a boson bath can be written in 
a general expression as 

/ 

P(t) = 



Poo(0) p 01 (0)e- B{At)+l£At/h 



(22) 



where e is the energy splitting between the electron energy levels. In short, pure dephasing 
cause a decay in the off-diagonal element of the density matrix for the two-level system that 
makes up the charge qubi t 21l36lS7 : 

Poi(t)~Poi(0)e- B2 W, (23) 
where the exponent function B 2 (t) is defined by 



B(t) = iki 1 d q ^f" sin 2 mth lSr- < 24 > 

Here u> q is the frequency of the phonons. In our study, we have investigated dephasing effects 
due to both acoustic and optical phonons. For acoustic phonons, we choose u; q = qc s for the 
relevant branches, while for longitudinal optical phonons, we choose cu q = wlo- It can be 
shown straightforwardly that the phonons that contribute significantly to pure dephasing 
are zone-center phonons (small |q| values), so that choosing linear and constant dispersion 
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for acoustic and optical phonons is an excellent approximation. The coupling constants g(q) 
due to deformation potential, piezoelectric and optical phonons are respectively given by 



#def(q) = D 
too(q) = Ml' (q) 



hq 



2pc s V 



X(q) 



n 



.x( q ) 



2pc s V 

M rr, \ 

^polar(q) = — ^=X(q) 



where X(q) is defined by 



X(q) = (r) 



^ + (r) 



(25) 

(26) 
(27) 



* + (r)\) = -^(q). 



(28) 



Here ± refer to the two states for the double dot charge qubit. Notice that all the integrals 
in this study are carried out using the Monte-Carlo technique. 



IV. RESULTS AND DISCUSSIONS 

A. Charge relaxation due to acoustic phonons 

We first calculate the electron relaxation rate from the excited charge qubit state (the 
first excited state of the double quantum dot) and explore its behavior as a function of the 
double dot parameters such as interdot distance and strength of the single dot confinement. 
Notice that throughout this paper, the QW width takes on a fixed value of 2L Z = 6 nm. 
We have also done calculations for a well width of 10 nm and the results are only slightly 
different. 

Figure Q shows the electron relaxation rate as a function of the interdot distance for an 
electron that is initially in the first excited state. The relaxation process is dominated by the 
emission of one acoustic phonon. For small inter-dot separations, the contribution due to 
deformation potential interaction is larger than due to piezoelectric interaction. For larger 
inter-dot separation, however, piezoelectric coupling becomes the dominant contributor be- 
cause of the different wavevector dependence in the deformation and piezoelectric matrix 
elements (^/q for deformation potential versus 1/y/q for piezoelectric interaction). Since 
phonon density of state goes to zero at small energy oc E 2 , the relaxation rate decreases 
with decreasing energy splitting between the initial and final states, which is the case for both 
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very shallow confinement wells and for largely separated quantum dots. Experimentally, it 
has been founcU^ that for large dots and large interdot distance (~ 150 nm) relaxation rates 
due to electron-phonon coupling should be smaller than 10 9 s~ l . Our calculations are quite 
consistent with these experimental observations. 




a (nm) 

FIG. 1: Relaxation rates of an electron in the first excited state through acoustic phonon emission 
as a function of the half interdot distance a. The relaxation rates due to deformation potential 
interaction, piezoelectric interaction, and the total relaxation rates are presented by dashed, dotted, 
and straight lines respectively. The strength of the lateral confinement is Tlujq = 3 meV. 

In Fig. El we plot the dependence of the relaxation rates on the confinement strength of 
one of the quantum dots (they are considered to be identical in this problem) for a fixed 
interdot distance of 2 a = 60 nm. The relaxation rates increase with increasing confinement 
strength as the energy splitting between the first excited state and the ground state increases 
(as illustrated by the inset of Fig. |3J), until the confinement strength reaches about hou ~ 
2.4 meV. Further increase in confinement strength causes increase of interdot barrier and a 
decrease in energy separation of the charge qubit states, so that relaxation rates also decrease 
as a consequence of the reducing energy splitting. The differences between the contributions 
of the two different types of electron-phonon interactions can again be interpreted in the 
same manner as in Fig. ^ The relaxation rates can also be given as a function of the energy 
splitting between the first excited state and the ground state (As) , as presented in Fig. 01 
Theoretically this graph more directly reveals the behavior of the relaxation rates: that it 
decreases monotonically with decreasing energy splitting between the initial and final states, 
basically because of the fast decreasing phonon density of state. In the double dot situation 
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FIG. 2: The relaxation rates versus the strength of the confinement. The scattering rates due to 
deformation potential, piezoelectric phonons and total relaxation rates are presented by dashed, 
dotted and straight line respectively. The interdot distance is 2a = 60 nm. 
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FIG. 3: Electron relaxation rates as a function of the energy splitting between the first excited 
state and the ground state in a double quantum dot. Again, rates due to deformation potential, 
piezoelectric interaction, and the total relaxation rates are represented by dashed, dotted, and 
straight line respectively. The energy splitting versus the confinement strength is given in the 
inset. The interdot distance is 2a = 60 nm. 

we study here for each energy splitting there could be two different dot configurations, as 
illustrated in the inset, thus there are two branches for Fig. El Notice that the energy 
splitting dependence of the relaxation rates is not universal,— because the electron-phonon 
matrix elements do sensitively depend on the form/size of the electron wavefucntions. 
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Since we have calculated the low energy spectrum of a horizontally coupled double quan- 
tum dot, we can easily calculate phonon emission rates when the electron is in an excited 
state. For example, Fig. 0] presents the phonon emission rates of an electron initially in the 
second excited state as a function of the half interdot dot distance a. Now a phonon-emitting 
transition can take the electron to either the first excited or the ground state. Furthermore, 
since the second excited state of a double dot is essentially made up of the 2p orbitals of 
the two single quantum dots, the energy splitting between the second excited and the first 
excited or ground states never goes to zero: at large interdot separation, this energy splitting 
approaches single electron excitation energy huo, which is chosen as 3 meV in this calcula- 
tion. Thus the phonon emission rate remains finite at large interdot separation, as clearly 
illustrated in Fig. In essence now electron relaxation is dominated by the electron-phonon 
coupling in each of the single quantum dots. This also explains why the relaxation rates 
into the first excited state and the ground state becomes identical as interdot separation 
increases. 



4xio" 




50 



a (nm) 

FIG. 4: The relaxation rates for the second excited state as a function of the half interdot distance 
a. The total scattering rates include contributions from both the deformation potential and the 
piezoelectric interaction. The straight line gives relaxation rates to the ground state, while the 
dashed line is the relaxation rates to the first excited state. The strength of the confinement is 
fuvo = 3 meV. 

In summary, our results on the electron relaxation through single phonon emission in a 
double quantum dot 38 show a relatively simple and straightforward dependence on energy 
separation of the initial and final states (basically related to phonon density of state), and 
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thus a sensitive dependence on the confinement strength and interdot distance (since energy 
splitting between the charge qubit states depends sensitively on these parameters). Accord- 
ing to Eq. El there is a strong angular dependence for the electron-phonon coupling matrix 
elements. However, integral over the phonon wavevectors essentially averages out all the 
detailed features in the final relaxation rates. Another feature of our results is that piezo- 
electric interaction dominates when the energy splitting between the charge qubit states is 
small, while deformation potential interaction dominates when the energy splitting is large. 
This feature again has a very simple physical explanation in the different q-dependence of 
the two types of interaction. 

B. Pure charge dephasing due to acoustic and optical phonons 

Electron-phonon interaction not only causes electron relaxation (or excitation at finite 
temperatures) between qubit states, it can also cause dephasing between them if the electron 
is in a superposition state, as we have discussed in Section III.C. Here we calculate the 
dephasing effects from both acoustic and optical phonons. According to Eqs. ([23]) and (JUJ), 
the quantity B 2 (t) completely determines the loss of coherence from the off-diagonal density 
matrix element between the charge qubit states, therefore it is the quantity we focus on in 
all the results and figures in the following. 

Figure El presents the time dependence of dephasing between the charge qubit states for 
an unbiased system due to acoustic phonons. According to Eq. PO]). dephasing should be 
quite small in the unbiased situation because it mostly comes from overlaps in between 
the double dot. An interesting feature of curves in Fig. El is that they rapidly increase for 
the first 10 ps or so, then more or less saturate, so that B 2 (t) depends only very slowly 
on time after 100 ps. Mathematically, the very fast time- dependence of dephasing is due 
to the trigonometric dependence on phonon frequencies and time sinu; s i/2. As indicated 
in Eq. (|20p . only zone-center phonons contribute significantly to dephasing (because of the 
integral over exp(iq • r)), with frequency ranging from zero up to ~ hc s /a, which is in the 
order of 0.1-1 THz for GaAs. As time evolves starting from zero, the zone-center phonon 
contributions quickly mix and the initial rise of B 2 (t) is mostly determined by the higher 
frequency phonons because their density of state is much higher. After the initial rise, 
phonons with different frequencies will not contribute always constructively, thus producing 
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the much flatter behavior of B 2 (t). Furthermore, B 2 (t) will not rise monotonically after the 
initial rise because the phonons contribute through a sinusoidal function. The time evolution 
shows that dephasing quickly rises but then saturates after about 100 ps. 




FIG. 5: Dephasing rates as a function of time t. The dashed, dotted, and straight lines represented 
dephasing rates due to deformation potential, piezoelectric interaction, and the total, respectively. 
The strength of the confinement Tiluq = 3 meV, the interdot distance is 2a = 60 nm, and there is 
no interdot bias: Vr = meV. 

The dephasing behavior here is quite different from the ordinarily assumed exp(— 7 P /ii) 
type of behavior. This difference is closely related to the spin-boson type of coupling in 
the present proble m 36 ! 37 and to a degree the phonon density of state of the semiconductor 
structure. There are two important consequences for the temporal behavior of B 2 (t). First, 
there is a very fast initial dephasing, occurring in a time period smaller than 100 ps, due to 
the interaction between the qubit electron and the acoustic phonon bath. Second, the time- 
dependence of B 2 at large time is very flat — it can basically be taken as a constant after 100 
ps. A constant dephasing factor will not produce a decaying signal in terms of, for example, 
oscillations in electrons. Instead, it simply reduces the contrast in the charge oscillation. 
This can be seen easily from Eq. 1221 The presence of a constant exp(— B 2 ) ~ exp(— 0.05) 
simply reduces the magnitude of poi by a constant factor of 0.05, which is not a particularly 
large suppression (though significant in terms of fault tolerant quantum computing). 

In Fig. El we further explore the dephasing rates as a function of the interdot distance for 
an unbiased system. As indicated in Eq. (ffilj) . if the two quantum dot are well separated, 
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the overlap integrals go down quickly, so that dephasing should also be strongly suppressed. 
This is exactly the behavior we observe in Fig. 




FIG. 6: Dephasing factor B 2 (t) as a function of the half inter dot distance a. No bias is applied 
across the two quantum dots, and time t is chosen to be 60 ps. The dephasing rates due to defor- 
mation potential, piezoelectric interaction, and total dephasing rates are represented by dashed, 
dotted, and straight line respectively. The strength of the confinement is Tlujq = 3 meV. 

When a bias is applied between the two quantum dots of the double dot, the magnitude of 
dephasing should increase, which is exactly what our numerical results presented in Figs. [7| 
and |H1 follow. Figure [7| shows a very similar temporal behavior as that in the unbiased cases, 
albeit with a much larger saturated value for B 2 . Figure |H1 shows the dependence on the 
interdot bias voltage Vr by the dephasing rate. As expected from Eq. (J20|) . the dephasing 
rates increases but then saturates as the qubit states at high bias are essentially the two 
single dot ground states (assuming higher excited states still have not influenced the charge 
qubit states yet). In Fig. [S] we have also included the contribution of optical phonons to 
dephasing. Although this contribution is smaller than the acoustic phonon effect, it is still 
a considerable contribution. 

As we have pointed out before, optical phonons do not contribute significantly to elec- 
tron relaxation. However, the zone-center longitudinal optical phonons do contribute to 
dephasing through the polar interaction. To better understand the pure dephasing due to 
the electron polar interaction with optical phonons, we can perform an analytical assessment 
for a simple double dot configuration. For optical phonons, u; q = ujlo is a constant near the 
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FIG. 7: Dephasing rates as a function of time t in the presence of interdot bias. Again, the 
dephasing rates due to deformation potential, piezoelectric interaction, and the total rates are 
represented by dashed, dotted, and straight line respectively. The strength of the confinement is 
huo = 3 meV, the interdot distance is 2a = 60 nm, and the interdot bias voltage is Vr = 1.5 meV. 

zone center. Thus at zero temperature 

M 2 . 2 iV LO t f , 3 |J(q)| 2 



B{t) ~ n^ LO y 8111 q q 2 



2e 2 / 1 1 \ . 2 ujLot 



L) sin 2 ^l f ^|J(q)| 2 (29) 

Cn / 2 J 



The g-integral is now just a number that is inversely proportional to the size of the double 
dot. The magnitude of B 2 is thus determined by the polar interaction strength and the 
double dot size, and the time-dependence of B 2 is all in the sinusoidal factor. Therefore, 
pure dephasing due to optical phonons sets in at a very small time scale, in the order of 
100 femtoseconds because of the fact that Hulo ~ 36 meV. The magnitude of B 2 is not 
vanishingly small, either. Using nominally GaAs parameters and assume spherical quantum 
dots with Gaussian wavefunctions ~ exp(— r 2 /2a 2 ), we can estimate the order of magnitude 
of the dephasing factor as 



^0.05 sin 2 (30) 

For the last step we assume a wavefunction size a ~ 20 nm, which corresponds to a pretty 
small quantum dot. Since B 2 is inversely proportional to the size a of the quantum dot 
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FIG. 8: Dephasing rates as a function of the interdot bias voltage Vr. The solid line represents 
dephasing rates due to electron-acoustic-phonon interaction through both deformation potential 
and piezoelectric interaction. The dashed and dotted lines represent the deformation potential and 
piezoelectric contributions separately. The dot-dashed line represents the dephasing effects from 
polar interaction with optical phonons. Here the strength of the confinement is Tiojq = 3 meV, the 
interdot distance is 2a = 60 nm, and the time of observation for the dephasing effect is t = 60 ps. 

wavefunction, larger quantum dot would produce a smaller dephasing magnitude for B 2 . 
Anyway, it is clear that optical phonons produce a dephasing effect that has a comparable 
if somewhat smaller magnitude as the acoustic phonons. 

The dephasing effect from optical phonons evolves extremely fast, so that the only ob- 
servable effect would be its average over time, which is a constant. This is quite similar 
to the pure dephasing effect from the acoustic phonons, which also rises rapidly (~ 100 ps, 
slower than optical phonons but still much faster than the ordinary time scale of nanosecond 
for charge dynamics). Recall that a constant dephasing factor will only reduce the contrast 
in the measurable quantities such as charge oscillation. For optical phonons this amounts 
to a reduction in the magnitude of p i by a constant factor of 0.95, which, like for acoustic 
phonons, is not a large suppression. Experimental techniques have not developed to the 
degree to allow detection of such a small suppression of signals. 

In short, dephasing effects on the electron orbital degrees of freedom from electro n-phonon 
interaction in a double dot should reveal itself mostly through a reduction of contrast in 
measurable physical quantities (such as electron oscillation between the double dot), but not 
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a temporally decaying signal. Decays observed in experiments such as [15| should originate 
from relaxation, not dephasing, if it is dominated by electron-phonon interaction. 



V. CONCLUSIONS 

In this study we have investigated electron decoherence in a double quantum dot due 
to electron-phonon coupling. In particular, we have evaluated electron relaxation through 
emission of a single acoustic phonon. We found a sensitive dependence of the relaxation rates 
on system parameters such as confinement strength and interdot distance. We have also 
evaluated electron dephasing through interaction with both acoustic and optical phonons — 
because of the absence of energy conservation requirement, all phonon modes contribute to 
dephasing. 
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